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Abstract 

The influence of the irregular forces on the evolution of a triple hierarchical stellar system moving in 
the field of stars have been studied. Triple hierarchical stellar systems are stable in contrary to the 
stellar systems with comparable distances between all components. We considered the motion in the 
frame of the general tree-body problem using differential equations of the motion with the Hamiltonian 
without short-periodic terms. 

For isolated triple stellar systems, where we took into account the perturbations until the third 
order, we obtained the solution in which the mean motions of both components have the secular 
accelerations. Under the influence of perturbations of the distant component the mean motion in 
the near pair is slowed and vice versa. The mean motion of the distant star is constantly increasing. 
These changes are small, but on the cosmological time interval the hierarchical systems will convert 
into stellar systems, in which all components have comparable distances between each other. Such 
systems become unstable. 

In a general case, if we take into account the irregular forces, the angular momentum of this system 
and its summary energy might be either loss or gain. These changes may influence on the dynamical 
evolution and stability of the stellar system. 

Introduction 

We have investigated the dynamical evolution of the triple stellar system using the analytical theory 
of the particular case of the nonrestricted three-body problem (Orlov and Solovaya, 1988). Masses of 
components are comparable and the ratio of the semi-major axes of their orbits is the small parameter. 
The value of the eccentricities and the inclinations of both orbits are without any limitations. Unlike 
multiples systems with the comparable distances between components, these systems are thought to 
be stable. 

The dynamical stability is understood as the conservation of the configuration of the system 
over the astronomical long time interval - the eccentricities of the orbits remain less 1, the mutual 
inclinations change in small limits, there are no close approaches among the bodies, and the triple 
system is located inside of a stellar cluster. The study of the isolated triple stellar system, where we 
took into account the perturbations until the third order in the Hamiltonian (Solovaya, 1977), showed 
that the mean motions of both components have the secular accelerations (Solovaya, 2010). 

Due to the perturbations of the distant component the mean motion of the close pair slows, while 
for the distant component increases. These changes are small, but on the cosmological time interval 
hierarchical systems will converted into stellar systems, in which all components have comparable 
distances between each other. 

If the triple stellar system is located inside of a cluster, it must be subjected to perturbations of 
the surrounding stars. We have done attempt to estimate the effect of their influence on the dynamical 
evolution and stability of the triple stellar system, moving through the gravitational field of surround- 
ing stars. In general case it may be either loss or gain of the kinetic energy. We used the Harrington 
S criterion of stability (Harrington, 1972). For direct motion, when the mutual inclination of the two 
orbital planes is less than 90°, S = a,2 (1 — e-i) / o-i > 3.5. For the retrograde motion S > 2.75. Here a\ 
and a,2 are the semi-major axes of the orbits and is the eccentricity of the distant component's orbit. 



The gravitational perturbations 



The motion of the isolated triple stellar system was studied using the analytical theory of the general 
three-body problem (Orlov and Solovaya, 1988). We have used the Hamiltonian without short-periodic 
terms. As the intermediate orbits we used nonkeplerian ellipses. For the computation was used the 
Jacobian coordinate system and the canonical Delaunay elements - Lj, Gi, k, gi (i=l,2) taking the 
invariable plane as a reference plane. Expanded in terms of the Legendre polynomials and truncated 
after the terms of the third order, the Hamiltonian has the form: 
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where R are the perturbations of the third order, the coefficients 71, 72, and 73 depend on masses, 
and 
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c is the constant of the angular momentum, g\ is the argument of the periastron of the close pair 
in the invariable plane, and q is the cosine of the mutual inclination of the orbits. The solution in 
hyperelliptic integrals 1% (i = 1, 2, 3) with the Hamiltonian until the terms of the second order was 
obtained by the method of Hamilton- Jacobi (Orlov and Solovaya, 1988). 

Denoting the mean motions of stars i/j, which equals 2 n P~ l , where Pi are the periods of revolution, 
we found that the mean anomalies l\ and £ 2 are expressed as: 



l\ = B\ + v\ (t — to) + periodic terms, 
h = B2 + V2 (t — to) + periodic terms . 
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where Ej (i = 1, 2, 3) are the secular parts of the hyperelliptic integrals. 

For stars moving along unperturbed orbits their mean motions may be expressed according to the 
third Keplerian law by formulae: 
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where ai and a 2 are the semi-major axes of their orbits. Due to the mutual perturbations the value 
of the mean motion of the close pair slows by Ani = v\ — ri\, while that of the distant component up 
by An 2 = v 2 — n 2 . 

To the illustration of the theory the stellar system e Hydrae (ADS6993) was selected. The three 
brighter components of the system are moving on the orbits with the semi-major axes a\ = 3.967 AU 



and a 2 = 75.76 AU. The eccentricities of their orbits have values e\ = 0.67 and e 2 = 0.29. The masses 
of the components in the solar mass are mo = 1.50, m\ = 1.30, and m-2 = 2.74 (Heintz, 1963). 
The results showed that for the close pair 

v x - m = -0.019° year" 1 , (7) 

and for the distant component 

v 2 - n 2 = +0.002° year -1 . (8) 



Adding the perturbations of the third order R in the Hamiltonian the mean motions v\ and v 2 
have the secular accelerations. The accelerations are small but on the cosmological time scale may 
change the configuration of the system. 

The perturbation function has the following form: 
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The mean motion v\ of the close pair changes with time according to the formula 

(3) 
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while the semi-major axis a± of the close pair increases. Mean motion v 2 of the distant star increased 
with time too, according to the formula 
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where o\ and a 2 are the secular accelerations. Finally, the semi-major axis a 2 of the distant component 
decreases. 

The evolution of the semi-major axes a\ and a 2 within the interval of 4.3 x 10 5 years under the 
influence of the gravitational perturbation is shown in Figure 1 and 2. It is seen that under the 
influence of gravitational perturbations the triple hierarchical stellar system converts into the stellar 
system in which the components have comparable distance between each other. The coefficient S is 
less than 3.5. Such system can be unstable. 



Effect of the dynamical friction 

Let the triple stellar system moves through the gravitational field of the surrounding stars in a stellar 
cluster. The triple system must be subjected to action of the irregular forces, due to random star 
approaches. One of them is the dynamical friction. In general case, the triple stellar system may loss 
or gain the kinetic energy. 

Stars in the cluster are considered as point masses interacting with each other according to the 
Newton's law. Each approach of the test star with the field star is based on the two-body approxi- 
mation. Mass of the test star is bigger than mass of the field star. We have no information about the 
motion of the stars in the system. We attempt roughly to estimate the influence of dynamical friction 
on the motion of the stars in the triple system located in an isolated globular or open cluster. Such 
clusters contains several thousand stars. We shall assume that the velocities of the stars in the triple 



system not differ too much from the mean velocity of the field stars. Then the effect of the dynamical 
friction may be expressed by the following simplified equation: 

Av 

— = —rjv , (13) 

where n is the coefficient of dynamical friction, Av is the mean change in the velocity of the test star 
during the time interval At. 

Chandrasekhar (1943) concluded that the coefficient of dynamical friction must be of the order 
of the reciprocal of the time of relaxation of the system. Consequently the effect of the dynamical 
friction can be apparent in systems with relatively short time of relaxation. We do not know the 
velocity of the every star. But we can use the virial theorem, which confirms that for the stable state 
in the gravitating spherical distribution of equal mass objects, the average potential energy must be 
equal the average kinetic energy, within a fractal two. 

From the virial theorem 

2<K tot >=<U>, < V 2 >=G n ™±, (14) 

tic 

where < K to t > is the average value of the kinetic energy and < U > is the average value of the 
potential energy. < V > is the average value of the velocity of a field star, R c is the radius of the 
cluster, n is number of stars in the cluster, and mj is mass of the field star. 

Introduce the dynamical friction relaxation time Tf r defined by the equation 

T fr = ^^T d , (15) 
rrif + nit 

(Bertin, 2000), where is the relaxation time relative to the deflection, defined as 

T < (16) 



and Td is the natural crossing time 



Then the dynamical friction equation is 

Av v t 4irG 2 m t p f F (v) In A 

-At=-f7 r = ^ V ^ (18) 

(Bertin, 2000), where nit is the mass of the test star, mj is the mass of the field star, F(v) is the 
relation of the velocity Vf of the field star to the velocity vt of the test star, pj = m^n is the mass 
density of the cluster. 

The Coulomb logarithm is defined as 

lnA = ln^, (19) 

where b\ and 62 are the minimum and maximum values of the impact parameters. 

We used formulae of the two-body approximation. If the velocity of a star increases (decreases) 
on the Av in its apocenter (pericenter), then the distance of the pericenter (apocenter) grows (falls) 
to the value 
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where p = k 2 (mo + mi). 



Results 



We have used Equation (18) for the estimation of the change of velocities of components in the triple 
stellar system in two cases, when this system is located in the globular cluster and in the open cluster. 

It is known that the average linear diameter R c of the globular cluster range from 20 pc to 2000 pc. 
For our calculations we adopted values n = 10 6 , R c = 35 pc, and G = 4.3 10 _3 pcMq 1 km 2 s -2 . For 
the globular cluster with the equal star masses m/ = 1 Mq we see result in Figures 1 and 2. In this 
case the dynamical friction have no the influence on the evolution of the semi-major axis of the stars 
in the triple system. 




Figure 1. The evolution of the semi-major axis of the triple stellar system in the globular cluster with 
n = 3 x 10 3 members during the interval t — 4.3 x 10 5 years under the influence of gravitational perturbations, 
dynamical friction and combined perturbations for the growing values of Av. 




Figure 2. The evolution of the semi-major axis of the triple stellar system in the globular cluster with 
n = 3 X 10 3 members during the interval t — 4.3 x 10 5 years under the influence of gravitational perturbations, 
dynamical friction and combined perturbations for the slowing down values of Av. 



For the next example we took the open cluster similar to the Pleiades, with the number of members 
n = 10 3 , the mass M = 800 M©, and the radius of the cluster R c = 2pc. But the impact parameter 
b\ in Coulomb logarithm was taken as a free parameter. The values of the impact parameters were 
selected b\ = 20 AU and 62 = 2pc. The gravitational constant G = 4.3 10 _3 pcMq 1 km 2 s -2 . 

The evolution of the semi-major axes are presented in Figures 3 and 4. It is seen that within the 
interval t = 5.6 x 10 5 years the coefficient S is greater than 3.5. If close approaches are possible, the 
dynamic friction has the influence on the evolution of the semi- major axis of the stars in the triple 
system. 




Figure 3. The evolution of the semi-major axis of the triple stellar system in the open cluster with n = 3 x 10 3 
members during the interval t = 5.6 x 10 5 years under the influence of the dynamical friction for the growing 

values of Av. 
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Figure 4. The evolution of the semi-major axis of the triple stellar system in the open cluster with n = 3x 10 3 
members during the interval t = 5.6 x 10 5 years under the influence of the dynamical friction for the slowing down 

values of Aw. 



Conclusions 

The mean motions of both components in the triple stellar system under the influence of the grav- 
itational perturbations have the secular accelerations. The mean motion of the close pair is slowed 
and vice versa. The mean motion of the distant star is constantly increasing. On the cosmological 
time scale the configuration of the triple stellar system converts from the hierarchical system into the 
stellar system, in which all components have comparable distances between each other. 

The influence of the dynamical friction on the change of the semi-major axes of the stars in the 
triple stellar system is more less in the comparison with the gravitational perturbations of third order. 
In the case of a close approach of the field star to the triple stellar system the strong perturbations 
are possible. In this case the perturbations from the dynamical friction are greater than the gravita- 
tional perturbations of the third order. The influence of the dynamical friction on the evolution of the 
semi-major axes of the both orbits in the triple stellar system is evidently and the system can stay 
hierarchical. 
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